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Abstract 

We study complex bosons and fermions coupled through a generalized Yukawa 
type coupling in the large-iVg limit following ideas of Rajeev[Int. Jour. Mod. 
Phys. A 9 (1994) 5583]. We study a linear approximation to this model. We 
show that in this approximation we do not have boson-antiboson and fermion- 
antifermion bound states occuring together. There is a possibility of having only 
fermion-antifermion bound states. We support this claim by finding distributional 
solutions with energies lower than the two mass treshold in the fermion sector. 
This has important implications from the point of view of scattering theory. We 
discuss some aspects of the scattering above the two mass treshold of boson pairs 
and fermion pairs. We also briefly present a gauged version of the same model and 
write down the linearized equations of motion. 



1 Introduction 

Quantum field theories are both fundamental and challenging. Despite the fact that our 
description of the world of elementary particles is based on quantum field theory, we still do 
not have a complete understanding of interacting field theories, especially their bound state 
structure. For this reason it is interesting to study simple examples where one can actually 
make more progress. There are various such models in two dimensions and they have been a 
valuable source for new ideas and testing ground for many years (see [0] for a comprehensive 
selection of topics). 
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In this article we study a two dimensional model which could be another possible toy 
model for understanding Yukawa coupled field theories in four dimensions. The physically 
important one is the gauged Yukawa theory, as we know from the present day version of 
the standart model. In this work we will also study some aspects of the gauge coupled 
interacting bosons and fermions. We will actually present the linearized equations of the 
gauged version at the end of our paper, but our main emphasis is to understand the simpler 
case without the gauge potentials (of course it is not so clear if this is really a simpler 
theory). A more complicated version of the model we discuss is investigated in |0] using 
path integral techniques. The results of the present paper are somewhat different since 
we follow a Hamiltonian approach (and we are not taking the most complicated possible 
version) . 

Let us comment on some fundamental work in the literature that we are aware of: the 
literature on Yukawa theory is vast, we will mention only a few of them, a rigorous construc- 
tion of two dimensional model is given in the papers 0, Q , the contstruction of the proability 
measure within the Euclidean formalism is done in |Q. The most recent rigorous analysis 
was given in by following a renormalization group type idea essentially inspired from [0. 
It will be interesting to attempt such a rigorous approach for the model we discuss below. 
The standart Yukawa coupling in two dimensions in the light-cone approach is discussed in 
P, 1^ . A further analysis using the Tamm-Dancoff approximation in the light-cone is pursued 
in 1 10 1 . A very interesting discussion of the Yukawa model in four dimensions is presented 
in [|11|] , clearly a four dimensional model has many more interesting features. A more recent 
analysis of the fermion bound states of the same model is discussed in [jT2[. The equiva- 
lence between the light-cone and covariant perturbation theory is analysed in |TB[. (Good 
reviews of field theories in the light cone are given in |]T^ and renormalization of light cone 
Hamiltonians in [|15|). 

In this work, we introduce a color degree of freedom for the purpose of reaching a Hartree- 
Fock type approximation. Following the ideas suggested by Rajeev in |]16| we reformulate 
the problem in terms of color invariant bilinears (see also Ijl^ for some similar ideas). The 
details of this refromulation are explained in our previous work within the context of 
gauge theories, therefore in this work we will use the results of the cited article directly. In 
some sense the present article is a natural continuation. The reader should consult ||18[ for 
more information on the geometry of the resulting classical phase space. We also recommend 
the lectures notes of Rajeev on two dimensional QCD [|19|. In [|16], it is shown that the 1-1-1 
dimensional QCD in the large- A^^^ limit can be reformulated as a classical field theory with an 
infinite dimensional phase space, which is identified to be the restricted Grassmanian. The 
study of two dimensional QCD in the large- A'c limit is given by 't Hooft in his well-known 
paper pO|, its generalization to scalar fields is done in ||21[ and in using Hamiltonian 



methods. The two dimensional combined fermions and bosons QCD is given by Aoki in 
2^ and also discussed by Cavicchi 0. Within Rajeev's approach one can reach the 



same results by using a linear approximation to the full theory. One can further study 
baryons by using a variational ansatz which does not correspond to small fiuctuations of the 
fields, therefore cannot be seen by the linear approximation. Here we obtain the general 
Hamitonian in the large-A'c limit for gauge coupled bosons and fermions which are also 
interacting through a generalized Yukawa type interaction. The meson equations are given 
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for the linearized theory. Our presentation is incomplete since we do not study beyond the 
linear approximation and we plan to come to a more detailed analysis in the future. 

We note that in the simpler model we discuss one can actually solve the integral equa- 
tions, ending up with some eigenvalue or scattering solutions. These equations require a 
simple renormalization to be meaningful (which at the end amounts to defining the singular 
integrals as the Hadamard principal value). We warn the reader that the form of the result- 
ing Hamiltonian suggests that a physically more relevant approximation in the non-gauged 
models could be given by a variational ansatz. This is due to the fact that the interactions 
in the linear approximation are all multiplied by the fermion mass. For heavy fermions we 
expect that the linear approximation gives valuable results, but for example in the case of 
massless fermions all the information is contained in the higher order terms, which cannot 
be accessed by the method we use. 



2 Coupling between Complex Bosons and Fermions 

We start with the action of our model with two Yukawa type couplings, 

(1) 

Here a refers to a common flavor index and it goes from 1 to Nc (we continue to write it 
as color symmetry, since at the end we will also talk about the gauged model). It is more 
natural to keep fiiY = fi2Y when there are no gauge fields, since we will use the color degrees 
as a way of reaching a mean-field description, but we will keep this more general form. We 
rewrite the action in the light-cone coordinate system, and we choose as our "time" 
coordinate, 

S = J dx+dx- (i^i)*j^d^t^La + iV2ilj*j^d+ipR^ + 0*"(-29_)9+0« - m^B0*"0„ 
There are many good introductions to the light-cone field theory, we refer the reader to 



25| , pq , p7| . We note that the lefthanded components of the fermion field are nondynamical, 
therefore we will remove ipia and its complex conjugate through the equations of motion in 
the quantum theory (we refer to [|I8| for our conventions in quantizing this theory, below we 
summarize the results), 

= 2id~^^^^ ' ^'^^^^^^^ ■ "^'"21- : (l>^'^(l>i3i'Ra ■ +mFi'Ra]- (2) 

We do not really need to worry about the normal ordering in the first factor, since in the 
large- A'^c limit these corrections will be of smaller order, we wrote it to emphasize that the 
reduction should be performed at the quantum level. At the second term we have a normal 
ordering for the color contracted bosons only. The second thing we notice from the light-cone 
action is that we are already in the Hamiltonian formalism. Therefore we can read off the 
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Hamiltonian directly from the action when we insert the solution of the nondynamical field 
back into the action. Thus we arrive at the following Hamiltonian, [from now on we write 
ip for ipji since this is the only fermionic field we have), 



+^f^lY ■ i^^i. ■ ^^"^V'a : 0^^0/3 : )• (3) 

This Hamiltonian as it stands is not normal ordered, to define it properly we need to normal 
order the color singlet products of bosons in the sixth term and the products of fermions in 
the last three terms. All these terms except one will give some divergences which can be 
cancelled by redefinitions of niBo and A^q in the original Hamiltonian as we will see. One 
of them cannot be removed by the original terms in the action and we will add a counter 
term which cannot be put into the original action. This Hamiltonian could be a better two 
dimensional representative of the four dimensional Yukawa theory, since in four dimensions 
phi-four coupling is necessary to renormalize the Yukawa interaction. 

Let us recall the quantization of this system in the light-cone coordinates, the Fourier 
mode expansions read. 

The normalization factors are chosen to reproduce the correct classical limits. To precisely 
define these expansions, we assume that the momenta range between (— oo, — eo] and [eo, oo), 
at the end of our calculations we set eo — > 0. This is physically meaningful due to charge 



conjugation invariance, and amounts to the principal value prescription, (see |jT8|, |T9[ for 
details) 

[X.b), = €5[V - q] , Kb), a^^m = sgn{p)6^J[p - q] . (4) 

One defines a Fock vacuum state |0) by conditions, 

aaipm = Xaipm = for p > at-(p)|0) = x^"(p)|0) = for p < 0. (5) 

(recall that we are assuming that there is an infinitesimal hole around p = to be taken to 
zero at the end of the calculations). The corresponding normal orderings are defined in [|18|. 
It is useful to keep in mind the normal ordering rules of the bilinears, 

: a-t(p)a^(g) := a"t(p)a^(g) - - sgn(p))5[p - q]. (6) 
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and 

: X^"(p)X/3(g) := X^"(p)X/3(g) - ^^1(1 - sgn(p))5[p - q]. (7) 

We provide some of the details of the reorganization of the Hamiltonian into normal 
ordered bilinears in the appendix. We formulate the theory in terms of the color invariant 
bilinears following the idea proposed by Rajeev |TB| and use our results in |T^. For the 



convenience of the reader we recollect some of the essential points: to define the large- iVc 
limit we introduce, 

M{p,q) = J- : : 

iV(p,g) = A : at"(p)a„(g) : (8) 

and their odd counterparts, 

q) = , Q{r, s) = ^a^-{r)xa{s) (9) 



In the large- Ai'c limit these bilinears become classical variables and we postulate 
the super Poisson brackets satisfied by these variables, which defines the kinematics of our 
theory: 

{M{p,q),M{r,s)} = -2i[M{p, s)S[q - r] - M{r,q)S[p - s] 

— 6[p — s]6[q — r]{sgn{p) — sgn{q))] 

{N{p,q), N{r, s)} = —2i[N{p, s)sgn{q)6[q — r] — N{r,q)sgn{p)S[p ~ s] 

— 6[q — r]6[p — s]{sgn{p) — sgn{q))] 

{Q{p,q),Q{r,s)}+ = -2i[M{p,s)sgn{q)6[q-r]+ N{r,q)6[p- s] 

+ S[p - s]6[q - r]{l - sgn{p)sgn{q))] 
{M(p,g),Q(r,s)} = -2t5[q - r]Q{p, s) 

{iV(p,g),Q(r, s)} = 2i5[p - s]sgn{p)Q{r, q) 
{M(p,g),Q(r,s)} = 2tS[p - s]Q{r, q) 

{7V(p,g),g(r,s)} = -2t6[q - r]sgniq)Q{p, s). (10) 

These classical variables now satisfy, 

M{p,q) = M*{q,p) N{p,q) = N*{q,p) Q{p,q) = Q*{q,p). (11) 

There are also constraints satisfied by these variables when we restrict ourselves to the 
subspace of the color invariant states. For our problem, this is an approximation, since there 
is no reason to expect that all the physical states are color singlets. In fact we will see that 
there are scattering states of our linearized equations. We write explicitly the constraints 
for the basic variables, 

{M + ef + QeQ^ = 1 
eQ^M + eQU + eNeQ^ + = 
MQ + eQ + QeN + Qe = 

{eN + eY + eQ^Q = 1. (12) 
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Above we use an operator notation, e(p, q) = —sgn{p)S[p—q], and {AB){p, q) = J dsA{p, s)B(s, q). 
The phase space of the resulting restricted theory is shown to be a super-Grassmannian in 
I8| , with its natural symplectic sturcture generalizing the results in [|1^ 



We can reexpress our Hamiltonian in terms of the above mentioned basic variables. After 
a somewhat long but straigthforward computation, the large- iVc limit Hamiltonian becomes. 



H = Hq + Hy, 



where. 



Ho = \mll^^Mip,p) + \ml^l 



\p\ 



N{p,p), 



(13) 
(14) 



[dpdqdsdt] 



5[p- q + s 



sq\ 



4- 



1 1 

p t 



+ K- 



Q{p,q)Q{s,t) 



^2Y^F f [dpdqdsdt] 



16 



6\p - q + s - t] 



rl 1 
- + - 
p q 



M{p,q)N{s,t) 



1 

64 



\pqkl\ 
[dpdqdsdtdkdl] 



\qskl\ 

IJ'lYfJ'2Y 



64 



6[p-q + s - t + k-l] 
1 1 



k — I — t s — t ~ q 



,,2 

64 V 



1 



+ 



iAy f [dpdqdsdtdkdl] 



64 



\pqkl\ 



6[p-q + k-l + s-t] 



s-t-l 
1 

t + l-k 



Qip,q)Qis,t)N{k,l) 
N{p,q)N{kJ)M{s,t).{15) 



Note that we have rescaled the coupling constants as fiiyNc ^— > /iiy,/X2yA^c ^ /^2y and 



X2 



As we discussed in the appendix, there are two possible renormalizations. In the first 
case we allow for the non-local counter terms and remove the divergent parts finding a local 
Hamiltonian-thus k, and X^b are just constants. If we only allow for the local counter terms 
then we find, 



A 



1 2 N , l^2Y_l^ 



s-t 



IT 



k-i 



l^R 



We introduce a renormalization scale ^n-, and assume that the renormalized values of the 
couplings, t^R^HR), X'^f^sil^R) '^ary with the scale fiR such that the Hamiltonian does not really 
depend on this scale. This means we should impose. 



kr{hr) = KR{iiR) - ^ In 

647r 



IJ'R 



l^R 
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2 



ALW = AL(/^ii) + ^ln ^. (16) 



TT 



The sign of the new coupling kr should not be fixed since it does not exist in the original 
action, and it has dimensions of mass square. 

For the rest of this work we will take the simpler Hamiltonian, that is we assume that all 
the renormalized couplings are ordinary numbers. Prom a more conservative point of view 
the local counter terms should be the general class of models we should investigate. We hope 
to return to a more detailed analysis in the future. 

This Hamiltonian along with the Poisson brackets and the contraint define our model 
completely. As it stands this is a complicated system. We plan to study a variational 
approach to this model in a future work. We will study a linear approximation to this model 
in the next section. 



3 The linear approximation 

We assume that all the basic variables deviate from the vacuum by small amounts, therefore 
we keep everything to first order. This means the linearization of the constraint and the 
linearization of the equations of motion. The constraint implies that 



M{u, v)=0 N{u, v)^0 Q{u, v)^Oifuv> 0. 
The equations of motion for u > 0,v < 0, found from 

dO{u, v) 



^{0{u,v),H}, 



(17) 



(18) 



where O refers to any one of our variables, could be linearized. Let us write these linearized 
equations of motion for all the variables. 



dM{u,v) 



.rup r 1 



1 



u V 



M{u,v) - i 



.ll2YmF 



1 1 
- + - 

U V 



N{s,s- (u-v)) 
s{{u — v) — s)\ 



[dsi (19) 



dN{k, I) 
dx+ 



and 



.m 



— i- 



BR 



1 1 
k~~l 



N(k,l)-i 



.IJ,2YfnF 



^J\kl 



[ds] 



- + 



s s — {k — I) 



8 "^J J\skl{{k-l)-s) 



..N{s,s-{k-l)), 



M{s,s-{k-l)) 

(20) 



dQ{u,v) 
dx+ 



RB 



U 



+ I 



Q{u, v) — K 

[dp] 



Ai f Q{p,p- [u - v)) 



\J\v{j> — {u — v) 
1 1- 



[dp] 



\{p — {u — v))v\ 



+ 

p u 



Q{P:P- (U-V)). 



(21) 
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The equation of motion for Q can be found by complex conjugation and does not carry 
new information. We note that the equations of motion for M and N are coupled, but the 
equations of motion for Q within the linear approximation is decoupled from the rest. 

So we will start with this one and make an ansatz as in |T^. Let us assume that 
the solution can be written as Q{u,v]x'^) = CQ{x)e^^+^~^ , where x = u/P^, P_ = u — v and 
define an invariant mass Aq = 2P_P+. (Strictly speaking we could take the solution of the 

form CQ{x)fQ{P^)e'^^+^~^, the arbitrary function fq factors out in the equations). Then we 
find. 



+ 



m. 



X 



+ 



m 



BR 



X 



4 Jo 



cq{x) - 8k 
dy 



dy 



^{l-y){l-x) 



y){l-x) 



In 



y 



)cQ{y)- 



This innocent looking equation actually requires a renormalization, as we will see shortly. 
Let 

-CQ{y) = A t^—CQ{y) = B. (22) 



I 

^{l-y) 



Then solve for cq{x)^ 



1-x) 



y^{l-y) 



{aA + hB)x + hA 



A^x(l 



X] 



m|(l 



x) 



m^x 



(23) 



where 



-8k, 



(24) 



A straigthforward solution will actually produce a divergence, the integration defining B is 
divergent. To find a finite result we need a renormalization prescription. Let us assume 
that the phase is given by e^iP++^P+i'^o))x+ ^ where 6P+{eo) denotes a divergent phase of the 
solution that we remove from the equations, and eo denotes a low momentum cut-off. The 
time derivative will drop a factor of 5P+(eo), and multiplying by 2P_ we denote it as SAq^eo) 
and rewrite the same equation as. 



m 



BR 



X 



1 



X 



aA 



hA 



X 



X 



+ 



X 



bBjeo) 



(25) 



Since the divergent part comes from the B term we expect that 6AQ{eo)cQ{x) can be taken 
as a counterm on the other side of the equality with the leading form — ^^M. The unknown 
function now is given by the same formula with a shifted coefficient of x, 



cq{x) 



x\ 



(aA + bB{eo) - ac{eo))x + bA 



Alx{l 



x] 



m. 



x] 



m^x 



(26) 
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Let us insert this back into (|22| ) and find the constants A,B{e()). After some algebra, we 
reach, 



A 



aA + bB{eo) — ac(eo) r/aA + bB{eo) — a^eo) 



2M 



(- 



2A2 



N bA , , ^ bA ^ 



m. 



2ml 



r bA 
2^ 



{Al-mls + ml) + bA\F,{AQ), 

(27) 

{Al-ml^ + ml) + {aA + bB{eo)-a,{eo))]F^{AQ), (28) 



where. 



F2 

fMq) 



Alx^ - {Al - m^B + ml)x + ml 



In 



dx 



lo A|x2 - (A| - mis + m|)x + mj, ' 
If we are looking for a bound state of a boson and a fermion, this requires, 

\mRB - mpl < Aq < muB + mp, 



then the last integral gives. 



FA An) = — ( arctan 
u ^ 



u 



arctan 



u 



where. 



^ = yi^Q- i^RB - mpy){{mRB + mFy - A^). 
Let us impose the two conditions, 

bA 

B^ = B{eo) 2" In(eo), bB{eo) - ac(eo) = bB^, 



m. 



then we see that if we set 



ac(eo) 



b^A 



In(eo), 



we can take eg — * 0^ limit, and keep B^:, A finite. The renormalized equations become, 



A 



aA + bB^ 



2Ai 



aA + bB^ 
2A2 



{Al - m^B + ml) + bA Fi(Aq), 



bA 
2ml 



{Al - m^B + ml) + (aA + 65,) Fi(Aq). 



(29) 
(30) 

(31) 
(32) 

(33) 

(34) 

(35) 
(36) 



If we solve for the ratio A/B^^ after some algebra, this gives us a consistency condition for 
the excitation energy Aq, 



F2-s(Aq)Fi(Aq) 



2aml 



Fi(Aq) F2 + s(Aq)Fi(Aq) 
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1 + 6Fi(Aq) F2 + {s{Aq) + ^)Fi(Aq) + 



2Af 



2im?pa 



(37) 



or equivalently, 



In 



m 



RB 



nil 



s(Aq)Fi(Aq) + ^Fi(Aq) 



In 



ml 



+ s{Kq)F,{Aq) 



1 + ^A*iymFFi(AQ) 



In 



m 



RB 



nil 



+ («(Aq) 



/xiymML . . A^ 



128k 



,(38) 



where s(Aq) = A^ 



m 



RB 



mi 



This equation is written in terms of dimensionless ratios of 



the variables, and it should be investigated numerically under the conditions we have stated 
before for Aq. Instead of numerically solving these equations we will look at one extreme 

case, when Aq niBR + nip (weak couphng). (It is interesting to investigate the opposite 
hmit Aq ^ \nip — m^^l, but the result is not so simple to interpret). It is better to use a 



different variable to study such limiting cases, we define A, via Aq 
Note that we have — 1 < A < 1. Our function Fi(Aq) becomes, 

1 1 / . r oj + A 1 



m 



BR 



+nip+2niBRniFA. 



Fi(A) 



mpniBR Vl - A^ 



arctan 



+ arctan 



VI- A2 



(39) 



Here uj = msR/niF and if we take A — > 1~, that means Aq — >• {niBR + nip)' and A — >• —1"*" 
corresponds to Aq \niF — niBRY' ■ There is nothing subtle about the first limit, keeping 
everything to first order gives us, 

-2 



1 - 



l-(2- 



K 



niF + niBR 



r 



nip + niBR 



In 



m 



BR 



W?p 



(40) 



We assumed all the way A fa 1, this could be consistent if for example we choose the couphng 
constant /xiy such that /xiy « niB when the ratio mp/niBR is not too large and if k < 0. 
There are other possibilities, but this simple one shows that there are solutions with the 
expected behaviour. If we assume A — > — 1+, the function Fi(A) — > l/mpniBR, and the 
calculations are more complex. If we set A = — 1+5^, we have Fi(— 1+5^) 



1 



This can be used to study the opposite limit, but due to its algebraic comlexity we leave it 
out, and only state that for various possible cases to be consistent, we find that nip » niBR 
is a necessary condition. 

Following the same strategy, we look for stationary solutions for the coupled equations: we 
start with the following ansatz for M, N: M{u, v) = CM{x)e^^+^^ and N{k, I) = CN{y)e^^+^^ , 
where x — u/Pm, Pm = u — v and similarly for N. If we now substitute these into the 
coupled equations of motion we see that the oscillations in time cancel out (since we select 
the same P+) and we end up with. 



nip 



1 



X 



X 



IJ'2Ymp 



X 



X 



m 



BR 



1 1 



X 



Cn{x) - 



fl2YmF 



An^{l-x)x Jo 



^y{l-y) 



I dy(- - -^)Cm 
Jo 1- 



dy, (41) 



(y) 
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+ 



where we set = 2PmP+, and = 2P/V-P+ for the invariant masses of the excitations 
of M and respectively. We notice that the desired decouphng of the total momentum 
variables would not have happened in the above equations if we had used the more general 
Hamiltonians. 

Before we plunge into the standart way to solve these equations we will talk about an 
interesting possibility, if we admit that the distributional solutions are also proper solutions 
of these equations. 

Let us write the equation for M, if we call 



then we can rewrite it as, 



Jo Jnih _ n,\ 



f^2Y 



^mCm{x) = rrip h kMKx), where d = . (44) 

L X 1 — X ^ 47:111 F 

For the solution we have in mind we need to impose 

dD < 1, (45) 

otherwise the energy will be unbounded from below. Let us assume that the last two terms 
in the equation cancel against each other. This condition implies that 



/' dy{^ - j^^Miy) = ff -^=Uyl (46) 
Jo 1- y^ 8tt Jo J„n _ ,A 



An Jo ""^y l-yJ^"'^"' Jo ^y[i-y)' 
or equivalently, 

^f'ly(--j^Yu(v) = ^fD. (47) 

47r JO ^y l — y^ oTT 

Thus we have to consistently choose everything to satisfy these conditions. The equation for 
M and A^ can be solved by using ^m{x) = 5{x — xp) and ^n{x) = S{x — xb)- What should 
we take as xp and xb? One way is to minimize the excitation energy for fermions, and then 
fix the bosonic parameter to have the cancellation. The eigenvalue for M becomes, after the 
minimizing choice is made. 



A^ = 



1 + dDy/^ + (1 - dDy/^] ' < Ami. (48) 



The last inequality is interesting since it implies that the fermions actually form a bound 
state at this energy. Now the condition we should have for the cancellation reads, 

l+dD^l/2 fl-dD\l/2 



l-dD^ ^1 + dD 
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where 



From here we can solve for D, 



f 



2lJ,2YmF 



1 



(50) 



(51) 



We should have > 0, this puts a condition on our couplings and the fermion mass. But 
there is a stronger condition, once we have the solution for the value of D, we can find the 
parameter to choose for the bosons from the definition of D, 



D 



dy 



-Uy) = 



1 



(52) 



It is possible to find xb H D > 2, since the minimum of the function on the right is 2, thus 
we need > 4. This implies a condition on our couplings, 



1 //^2y\ 



-V2 ^ A|b 



m 



2 ' 



(53) 



where we used dimensionless variables to express this inequality. 

We have the other condition about D, which says dD < 1. This is actually satisfied by 
our solution, so we need 

2<D<l/d. (54) 
From these we have a condition on the strength of the Yukawa coupling constant. 



//2y < 27rmF- 



(55) 



Now we can go back and find the actual value of the fermion bound state and the mass 
of the boson pair. The boson pair mass is simply given by 



m 



RB 



m 



RB 



xb{1 - xb) 
Similarly we have for the fermion pair. 



1 

d2 



4 
P 



16m 



RB 



Tt ni l 



Alt 



>4m^5. (56) 



rrii 



1/2- 



We note that the above solution is quite interesting: we assume that the relative strength 
of the Yukawa coupling is small, i.e. 112Y « Ar, and expand the square roots. 



2(i2 

A^«2m^(l + — ) 



2mi 



1 + 



1 ^^^Y 



^ A^ 



(57) 



This is not the result one should expect from a perturbative point of view. 
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Furthermore the boson pair becomes in this approximation, 



16 



TT 



-2 



m-BR > 16m^s. 



(58) 



If we further assume that 112Y « mp this imphes that the boson pair mass becomes very 
large. 

Actually, there is a whole range of solutions with iuix) = 6{x—xf) and ^n{x) = 5{x—xb)- 
We are free to choose one of them, say xp then the other one will be determined by the same 
consistency relation as above. We see that the boson pair excitation will always be bigger 
than 2'mBR, since the minimum is given by xb = 1/2. For the fermion pair we choose the 
consistent xp^s such that the mass is less than the two mass threshold. Let us briefly present 
our findings using the same notation as above. 

Let us search for a solution of the equation 



1-dD 1 + dD 



+ 



X 



X 



inix) = K]^^m{x) with A^ < Am]p. 



(59) 



We assume again that dD < 1, then the equality is satisfied if we set ^m{x) — 5{x — xp)., 
xp = 1/2(1 — adD) for < q; < 1. We require the same delta function solution for ^^v, this 
means we should solve for the equation, 

^ ^ (60) 



If we solve for D now, we find 



Xp 



1 — Xp 



Aad 



f 



1 



a 



(61) 



This implies that the first factor in the big parathesis should be positive. This is true if 
A|j > ^aiily Again to have a solution for xb we need 4 < D'^, this means 

Aad , , 

1 > 0, (62) 



Ao?d^ + 



/ 



which can be satisfied if the quadratic eqaution for d has real roots and we choose d in 
between-assuming alpha is chosen. This implies for \r an inequality. 



2a //i2y\2r _ / a\i2Y 



TT ^VTLp 



^2'Kmp 



< 



3. 

rup' 



(63) 



This means that we should choose Ar above a certain value, and this condition is stronger 
than the first one we found above. A uniform bound for various values of a can be chosen. 



2 /■fJ'2Y' 



27rm/ 



< 



3l 

ml 



(64) 



Incidentally this requires /i2y < 27rmp. We still have dD < 1 to satisfy. If we simply use 
this condition assuming a as given we arrive at the positivity of a quadratic expression in a: 

2a jily 



o? + 



TT \\ 



1 > 0. 



(65) 
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Notice that the range of allowed a will be bigger if we take [i2Y I^r ratio as large as possible. 
If wc use the uniform lower bound for A/j we will find the largest region, and if we denote 
the deviation from this value by a multiplicative factor A; > 1, we can insert this ratio into 
the quadratic expression, 



k 



a 



1 > 0. 



(66) 



The positivity is guaranteed if we choose a outside of the region between the two roots. The 
choice consistent with < a < 1 gives us. 




1- 



+ 4A;2 - 



1 - 



< a < 1. 



(67) 



The lower bound is a decreasing function of k, so the stronger relative values of \r will 
have a smaller domain of a's. Since for these choices we have a continuous range of a's, the 
spectrum of the problem is rather different. Such distributional solutions are not eigenvalues 
but they typically refer to the continuous part of the spectrum. In some sense these are still 
scattering states. This means we cannot use the free parts of the original Hamiltonian to 
study the scattering theory below the two mass thresholds. Above these values we will see 
that the scattering theory can be studied by conventional methods. Perhaps below this we 
need to use the minimum value of the spectrum to define new effective pair mass for the 
fermionic sector. We are not able to resolve this issue at the moment. 

If we now go back to the standart approach, again as in the case of boson-fermion 
pair, we will need to renormalize our equations, by assuming a divergent common phase 
5P+(eo), 2i^A^e*^^'^+*-'^°''^^+''^^, and similarly for N. The derivative will bring terms of the 
form 5AAr(eo) = 2P]\f6P+{eo) and 5AM(eo) = 2PM<5-P+(eo)- An inspection of the resulting 
equations show that, we have the leading behaviour 



SAM{eo)CM{x) 



Qc(eo)(l - 2a:) 
x{l — x) 



x(l 



(68) 



X) 



In general, /5c(eo) = cradeo), as we will see the precise value of a is not important. Now we 
can solve the unknown functions. 



;i - 2x){aA{eo) - adeo)) 
A\^x{l ~ x) — m'j 



^m{x) - V^(l ^) A%x(l-x)-ml. ' 



where 



and 



A{eo) 



dy / ^ 5(eo) 
^0 ^Jy[l-y) 



''BR 

60 y 



An Stt' 



(69) 
(70) 
(71) 
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If we introduce 



F(A,m) =V 

Jo 



dy 



Iq A'^y{l — y) — m?'' 
we find from the defining conditions of A(eo), -B(eo); 



(72) 



5(eo) = 
A{e,) = 

We now define aA^ 
into our equations, 



1 

aA{eo) - ac(eo)) — Ineo + (aA(eo) - ttc(eo))(^ - 2)F(AM,m^) 

777; 171 ^?72 TP 



2ml 



2aB{eo) + M(eo) - ^«c(eo)J i"(A^, ^^b)- (73) 
aA{eQ) — adeo), and -B,,, = -B(eo) Ineo; and insert these back 



If we set 



A2 



2)F(AM,m^) 



2aB^ + bA^ + —A^ In eo + [- 



mi 



aF{AN,mnB] 



]a,{eo) F{AN,mRB)- (74) 



ac(eo) = - 



a 

a 



-^A^ In eo, 



(75) 



we can take eo 
equations. 



aF{AN,mRB] 

0+ limit while keeping A^,,B^, finite. These will be our renormalized 



2ml 



- 2) F(Am, mp) A, = [2aB, + 6A,]F(A^, m^ij) 



(76) 



Incidentally we note that this physical prescription implies that the proper way we should 
define these integral equations is to use the Hadamard finite value (see the similar issue in 
TP| and 1^ for a recent discussion of the renormalization and distribution theory). 

Let us assume that we are looking for bound state solutions, the principal value integral 
in F{A, m) then becomes an ordinary integral. Now we have two different expression for 
the ratio B^/A^^ which give us the desired eigenvalues when we require square integrable 
solutions. If we assume that both of the eigenvalues are bound states we find. 



i4y V^m 



iml 



A 



arctan 



M 



A 



M 



4ml 



A2 



Stt 



An si Am 



''BR 



A2 



arctan 



A 



N 



A2 



-1. 



(77) 

Since the left hand-side is positive and the right one is negative this has no solution! if 
we demand both boson-boson and fermion-fermion pairs to form bound states, there is no 
solution. But we should not be alarmed, what we can do is to demand a resonance for the 
bosonic sector, then. 



Al^ - Ami 



47^2 



A 



arctan 



M 



A 



M 



Ami - Ai/ 



Stt 



A 



N 



Al - Ami J, 



In 



An + 



A 



N 



A^ 



1, 



(7J 
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and this has solutions in general. The resonance case seems to be special to 1+1 dimensions, 
the other possibility is to require bosons to have scattering states and fermions to be bound. 
This can be studied by analyzing the pole structure of the analytic continuation of the 
scattering amplitudes that will be worked out below. Due to the algebraic complexity of the 
resulting formulae, we will not be able to answer it in this work. 

We study the scattering states when we are beyond the two mass threshold. Let us make 
a digression for the moment and study a simpler problem, the lamb da-phi- four coupling. It 
is easy to see from our equtions for N{u, v) that the same ansatz for the solution leads to 

We assume that the operator on the right is acting on L'^{[0, 1]), with vanishing at the end 
points boundary conditions. The free part, 

2 

is an unbounded operator with a continuous spectrum [4m^, oo). This is easy to understand 
by studying a particle and antiparticle in the center of momentum frame. 

If the added term is not a "too strong" perturbation then the absolutely continuous part 
of the spectrum of the full operator on the right is the same as the spectrum of the free part 
and we can study the scattering states using the free part (there are various conditions we can 
state so that "too strong" becomes a precise statement, we recommend [0 for a thorough 



mathematical discussion of these issues, and |^ with more physical emphasis). The type of 



problem we study is analyzed in a recent valuable book by Albeverio and Kurasov . The 
interaction term is called a rank one perturbation. The Hamiltonian, 



m2 . A| 1 /•! 



H = —, - + -^^^^=/ dy , (81) 

X{1 - X) 87T l^^l _ Jo /^(l _ y) 



where everything acts on functions in i^^([0, 1]), can be written as 

H = H, + ^\f><fl (82) 
87r 

with < x\f >= f{x) = , ^ If such a perturbation is relatively form bounded then the 

X{\ — X) 

scattering states are given by the scattering states of the free part. To verify this condition 
it is enough to show that the added term satisfies 

^ T7l/ll<oo, (83) 



"(|i^o| + l)i/ 

where 1 1 . 1 1 denotes the usual norm. It is now simple to check that 



i< / dx\— -) ( , ) < 00. 84 

^ Jo yx{\-x)) ^ ' 
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In fact in the above problem we can find the resolvent of our integral operator (see |]32| 
{RH{Z)f){x) = [{H-Z)-'f]{x) 



(: 



\ -1 Ad ~ \ X\L — X f 

z) f{x)+^Aiz) y/ \ / 

^ QTT Zx[l — x) — Jo 



x{l — x) ' Stt Zx{1 — x) — m? Jo Zy(l — y) — m? 

for Z outside of the spectrum (and complex in general), and here we use the analytic con- 
tinuation of ^(A) to complex numbers and its explicit form is given below. The knowledge 
of the resolvent gives everything about the operator, for example we can find the spectral 
density function p(A) (which heuristically corresponds to the "eigenf unction" expansion) by 
using the well-known Stone's identity, 

p(A) = — lim ( i i V (85) 



In the scattering theory of finite rank perturbations has been worked out by using 
rigorous methods. We will only content with the result that the scattering theory makes 
sense beyond the bound state thresholds for both particles, and we can find the resolvents 
explicitly. For simphcity of our presentation we study the scattering theory by the standart 
methods in physics, and only find the wave operators. One can see that for a given value of 
= j^Y^ have two roots. 



1 m 
4 ~ A 



2. 



1/2 



(86) 



In physics we typically think of two particles approaching and then scattering off to infinity. 
An inspection of the kinematics of a particle and an antiparticle pair in the center of momen- 
tum frame reveals that A+ corresponds to the particle moving in the positive x^ direction 
(which we may take as "incoming" states), and A_ corresponds to the particle moving in the 
negative direction. From a physical point of view, the scattering data should give us the 
information about transmission and reflection of the pair. An equivalent description would 
be to find the wave operator, Q which maps (in general the projection to the absolutely 
continuous part of the spectrum of the original operator) the Hilbert space to the scatter- 
ing states (if we take the absolutely continuous part of the spectrum and use the spectral 
projections corresponding to these values) of the interacting Hamiltonian: 

e = ^f- (87) 

(In physics one typically uses which takes the wave functions at time zero and evolve 
them to positive infinity, this requires the +ie prescriptions in the integrals, we will see that 
for our problem it is more suitable to define the principal value one, this is why we use Q). 

To find the scattering amplitudes we rewrite the eigenvalue equation in Lippmann-Schwinger 
form, 

= - A) + ^ [A^ - -'^^ /' dy^dy^. (88) 



ill 



xil — x) 



^x{l-x) ^y{l-y) 
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We formally use the eigenvalue Sx{x) — S{x — X) of the free Hamiltonian. We can now solve 
for A(X) = -4^, 



1 + 



X2 



:ln 



A - VA^ - 4m2 



-1 



(89) 



where A^ = ^ 
Hamiltonian can be written as 



47rA^(A2 - Arn?) A + VA^ - Am? J ^A(l - A) ' 
. Thus the wave operator acting on the (formal) eigenfunctions of the free 



m 



m 



My) 



y{l-y) x{l-x)i ^x{l- x)y{l-y). 



\S{y-X), 



with A{X) = ^A(l - A)A(A). The left side of the expression gives us the distributional kernel 
of the wave operator. (If we are interested in incoing pair we could restrict ourselves to A+ 
values). We can expand an arbitrary vector into a series of the form / = dXf{X)dx{x), 
and we find 



dxn{x,x)f{X), 



(90) 



and this makes sense in general. Note that this result is exact (within the linearized large- A''^ 
limit) and we have a complete charaterization of the set of scattering states once A{X) is 
given. 

We will study the scattering states of the coupled equations beyond the bound state 
tresholds, that is when the energies are larger than both 2mF and 2mBR. The discussion 
preceeding these indicates that the free part and the interating coupled equations may not 
have the same scattering states (see the distributional solution we present). We follow the 
same idea as in the above problem and restrict ourselves to the heuristic Lippmann-Schwinger 
type approach. One can also find the resolvent exactly and verify the formulae below by 
more careful analysis. 

We solve for the scattering amplitude by using a renormalized Lipmann-Schwinger equa- 
tion: thus we have for the scattering, 

Cm{x;A) = S{x - Am) + P ^,^^^ - m^ ^^ ~ 2x)aA.{A) 

S{x - Xn) + ^T2-Ti — -1 —^x{l-x){aB,{A) + bA,{A)), 



^n{x;A) 
where A* (A), B^{A) satisfy, 
A (A) = 



A%x{l -x)- mlji 



+ {aB,{A) + A(A))F(A, niBR) 



Aiv(l — Xn) 

1 1 f A^ \ 

a(A) = — — + a^,(A) — -4 F(A,mf), 

Am -1- — Am V^f / 



18 



and F{A,'m) is the same function as before and we choose Max(4m|^, 4m|,) < A^, with 

2 2 

^ Am(T-Am) ~ Ajvli-Ajv) " ^'^^ ^'^ ^'^^^ parameter and call it simply X. If 

we solve for the scattering amphtudes, 



A (A) 



l-( 



+ X^V(l-4A(l-A))ln 



Vl-4A(1-A) 



X 



A(l-A) 



In 



1 + ^1-4A(1- A) 
v/l-4(mWm^)A(l-A) 



1 -1 



X 



^1 - ^mlJml)X{l - A) l + ^l-4(mW^|)A(l-A) 

1 



1-2A 



In 



^1 - A{mlJml)X{l - A) 



^Bii ^A(l - A) 47rm^ ^1 - A{mlJml)X{l - A) 1 + ^1 - 4(mW^F)A(l - A) 



Note that this result is written in terms of dimensionless variables. We can read off -B*(A) 
as well, 



B^X) 



1 - 2A 



A(l - A) 2TTmF 



1-4A(1- A)ln 



1- J1-4A(1-A) 



1 + ^1 - 4A(1 - A) 



A^X). (91) 



The reader can check that the above results actually reduce to the phi-four theory results 
we have found if we set fi2Y = 0. 

For the sake of completeness we will also present the scattering solutions for cq{x) vari- 
ables. Below, we use the same shorthand symbols A{X) and -B*(A) as in the bound state 
equation for cq: The renormalized scattering equations become. 



cq{x;X)^S{x-X)-V 



we should set 



A^x^ — (A^ — m%ji + m'p)x + rrip 



2 2 
^2 _ _|_ f^BR 



^/T^[{aA{X)+bB,{X))x+bA{X)], 

(92) 

(93) 



A ' 1-A' 

not surprisingly mp + rriBR < A < oo, and ^(A), B^{X) found from their definitions, when 
we simplify the result it becomes. 



A{X) = 



X 



X 32A2 



64A2 64A2 

iiiyrriF 



(A^ - + ml)F, + t^lI^F,) (l + i^lf^F,)" V, 



1 -1 



B^X) 
where 



(1 + 



(1 + 

1 



Fi) 'f2 



XVT^ 



8m F 



(A2-mL + m|)F, + |^F2)A(A) 



Fo^ln 



^BR 


, Fi = -ln 


rrip 


u 



(94) 
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here m = y (A^ — (mp — mBRY){I^^ — {mp + rnBRY). 

These define the wave operators of our model, as discussed in the simpler model of phi- 
four coupling. We see that the results are fairly complex expressions. One should study 
various approximate forms of these equations, and a numerical investigation of the poles of 
the amplitudes should give information about the bound states. 

Our present approach has one weakness, our results are nontrivial since the fermions 
have nonzero mass. This makes it sensitive to the sign of the coupling-we treated them 
as positive, but the results are valid if we simply assume them to be negative. Physically 
more interesting case would be to study the massless fermions. If we set m^? = all the 
interesting information we have is lost in the linear approximation, and we should go beyond 
this. This observation suggests that it is necessary to study some kind of variational approach 
to understand the system better. We plan to investigate this in the future. 



4 Gauged model 



We will assume in this part that the model is gauged by introducing SU{N^ Lie algebra 
valued gauge potentials and refer for our conventions to [|l8l. In the light-cone coordinates, 
and setting = 0, the action becomes. 



i.*a\2 



The restriction to the color invariant states in the gauge theory is actually necessary to make 
the Hamiltonian finite. For this theory the large- A^^^ limit should be a better approximation. 
Furthermore one expects baryons in this theory, the geometry of the large- A^^^ phase space 
should be useful to find a variational ansatz (see for a nice discussion of these ideas). 
Following the same reduction process, the Hamiltonian becomes 



H = Hq + Hy + Ht 



G, 



where, 



\p\ 



N{p,p), 



(96) 
(97) 



Hy is as given in equation (|T^), and the gauge contribution is exactly given in |T8 



G 



— J [dpdqdsdt] 



9 f 

+ — [dpdqdsdt] 
64 J 



ip - t)2 
1 

+ 



6[p + s-t- q]M{p, q)M{s, t) 



+ 



[dpdqdsdt] 



{p - t)2 {q - sf 
q + s 6[p + s —t — q] 



S[p + S-t- q]'-^±^±£±^Nip, q)Nis, t) 



\pqst\ 



Q{p,q)Q{s,t), 
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Above we rescaled our coTipling constants by a factor of as before and g^N^ i— > g^. Let 
us use exactly the same substitutions as before for the basic variables we have, and simplify 
the resulting equations into, 



Jo 



m[x) - — 



dy 



IJ,2YmF 



1-x 



^Niy) 



2/(1-2/) 



TT Jo {y — x)'^ 
--dy. 



and 



( 2 , 1 \^ I . 9^ dy {x + y){2-x-y) 



Att^J (1 — x)x 



Jo ^y l-y^ Stt Jo ^x{l - x)y{l - y) 



{y) 



(98) 



Again we see that the equations for and ^jv are coupled and they should be solved 
together. 



TT 



1 



1 — X 



dy 



2 — X — y 



- 8k 



dy 



^(l-y)(l-x) 



1 /■! 



(2/-^)V(i-^)(i-2/) 

dy 



cqiy) 



^{l-y){l-x)^^ 



(99) 



These singular integral equations can perhaps be investigated numerically. The linear ap- 
proximation could be a better one for the gauged model, since the effect of the gauge inter- 
action is to bring a singular operator. The non-gauged models require renormalizations, it is 
possible that the above equations will behave better due to the singular operators in them. 
Another important application is to find a variational ansatz for the baryonic solutions and 
have a linear expansion around these solutions. We plan to study these issues in more depth 
in the future. 
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6 Appendix: Reduction of the Hamiltonian 



In this appendix we will give some of the details of the reduction of the Hamiltonian to the 
desired color invariant products. Let us recall that the Hamiltonian is given in equation 
number (|15|). Let us start with the term, 



lO- 

When we write this in terms of the Fourier mode expansions, it becomes 
Ii\y f [dkdldtdpdqds] S[k — I + t — s + q — p] 



(100) 



\stkq\ 



t-k + l 



16 J J\qk\ 



sgn(s) — 1 



[ds] 



s{k — / — s) 



Notice that the divergent integral is isolated and a principal value regularization calculation 
shows that 



V / [ds] 



sgn[s) 



1 



In 



k-l 



where is an infrared cut-off. If we are only allowed to introduce local counter terms in 
the original action we should introduce a momentum scale fiR so that we can separate the 
momentum dependent part and purely divergent part of these type expressions: (this point 
is somewhat difficult to decide for this particular term since it is not possible to write such 
a term in the original action) 



V / [ds] 



sgn(sj 



1 



1 



s{k — / — s) 7r(A; — /) 



In 



k-l 



1 



7T{k-l) 



In 



k-l 



+ In 



(102) 



If we remove the divergent part by a counter term of the form 



IGtt 



In 



[dpdqdkdl] 1 



the finite term comes out to be 
[dkdldqds 



8Kij(/i 



R) 



In 



57r 



qk 
k-l 



k-l 



5[p-q + k- l]x^''{v)aMW^{k)xp{l). (103) 



^J'R 



q 



-X^^{k)a^{l)a^^{q)Xp{s)5[k -l + q-s]. (104) 



If we require the theory not to have a dependence on the arbitrary scale we introduced, it is 
natural to demand that the residual coupling to vary under a change of scale according to 



K-Ril^R) = f^Ril^R) - ^ In 



I^R 



fJ'R 



(105) 



The reader may be alarmed by the nonlocal expression in the interaction, but if we actually go 
back to the position space, the inverse Fourier transform gives a term, upto some constants. 



dxdy(asgn{x -y)-\n \hr{x - y)\sgn{x - y))i}'^°'{x)ij)a{y)4>^^{y)i)p{x), 



(106) 
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which has a logarithmic correction to the sign function. (This behaves worse for the short- 
distance than the coulomb potential |a; — but we should interpret sgn(a; — y) = if a; = |/, 
so there is no singularity at the short distance). In the text we will only consider the cases 
where these non-local terms dropped, or removed by taking them as part of the counter 
terms in the action. 

(Notice that this is the interaction one finds if we use a parity broken model as in | |18| |. 
The sign of the remaining interaction term is not determined since it is not in the original 
action it should be left as an arbitrary parameter). 

As another example we will discuss the term, 

^mFfi2Y J dx~ ^^""-^T^a ■ (p" <P<t (107) 

a Fourier expansion and removing the vacuum expectation value gives us a term 
f^2YfnF f [dpdqdsdt] 6[p — q + s — t] 



st\ s-t-q 



■ X^"(p)Xa(g) :: a^^{s)ap{t) : 
^^^I^ : a^^s)a,{s) : j[dp]^-^^, 

V [[d/-^^ = -2\n^, (108) 

where we have e/j and Ar as the infrared and ultraviolet cut-offs respectively. If we introduce 
a boson mass counterm of the form 



and 



: a^f'{s)ap{s) :, (109) 



4 ^ en J \s\ ''^ 

this term will be cancelled. 

The other terms are also done in the same way, and the rest is to collect all terms to find 
the Hamiltonian in terms of large- A^^c bilinears. 
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